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Abstract. We consider an energy conserving linear dynamics that we perturb by a Glauber dynamics with 
random site dependent intensity. We prove hydrodynamic limits for this non-reversible system in random 
media. The diffusion coefficient turns out to depend on the random field only by its statistics. The diffusion 
coefficient defined through the Green-Kubo formula is also studied and its convergence to some homogenized 
diffusion coefficient is proved. 

On considere un systeme d'equations differentielles lineaires couplees conservant une certaine energie et 
Ton perturbe ce systeme par une dynamique de type Glauber dont I'intensite varie aleatirement site par 
site. Nous prouvons les limites hydrodyanmiques pour ce systeme non reversible en milieu aleatoire. Le 
coefficient de diffusion depend de I'alea uniquement par sa loi. Nous etudions aussi le coefficient de diffusion 
dcfini par la formule de Green-Kubo et montrons la convergence de celle-ci vers un coefficient de diffusion 
homogeneise. 

1. Introduction 

The derivation of hydrodynamic hmits for interacting particle diffusive systems in random environment 
has attracted a lot of interest in the last decade. One of the first paper to consider such question is probably 
[lOj where hydrodynamic behavior of a one-dimensional Ginzburg-Landau model in the presence of random 
conductivities is studied. In |19) . a lattice gas with random rates is considered and a complete proof of 
hydrodynamic limits has been given in [6], [20) . Other systems have been investigated such as exclusion 
processes and zero-range processes ([II [8j [U [121 El US])- Interacting particle systems evolving in random 
media are in general of non- gradient. Roughly speaking the gradient condition means that the microscopic 
current associated to the conserved quantity is already of gradient form. Otherwise the general non-gradient 
techniques f|15), |22| ) consists in establishing a microscopic fluctuation-dissipation equation which permits 
to replace the current by a gradient plus a fluctuation term. But, if the system evolves in a random medium, 
such a decomposition does not hold microscopically because the fluctuations induced by the random medium 
are too large, and it is only in a mesoscopic scale that this fluctuation-dissipation equation makes sense ([6], 

m)- 

In |121 114) , by extending some ideas of |16| , a simpler approach is proposed. The idea is to introduce a 
functional transformation of the empirical measure, which turns the system into a gradient-model, in such a 
way that the transformed empirical measure is very close to the original empirical measure. The advantage 
of the method is that it avoids the heavy machinery of the non-gradient tools but is unfortunately restricted 
to specific models. Even if the techniques developed in [B], [50] seem to be more robust than the precedent 
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approach, it is not clear that in some situations, as in the situation considered here, they can be apphed 
without a substantial modification. 

The interacting particle system we consider is the following. To a simple energy conserving linear dynam- 
ics, flips with site dependent rates are superposed. Fix a sequence {'^x)x of positive numbers and denote by 
{il{t))t>o the Markov process with state space and generator given by 

(1) {CfM^{Af){v) + iSf)iv), /:M^^M 

where 

iAf){r]) = ^(?7x+i - Vx-i)dr^J 

and 

with 77^ the configuration obtained from rj by flipping rjx'. {il^)z = "Hz '^^ z ^ x, {r]^)x = —"Hx- This system 
conserves the energy ^^Sx-, ex — ilil'^) and the product of centered Gaussian probability measures with 
variance T > are invariant for the dynamics. 

Let (7x)a; be a sequence satisfying ^ and PT|) . For example, the sequence {^x)x is a realization of i.i.d. 
positive bounded below and above random variables with positive finite mean. We show (cf. Theorem [1]) 
that, starting from a local equilibrium state with temperature profile Tq = l//?Oj the system evolves in a 
diffusive time scale following a temperature profile T, which is a solution of the heat equation 



(2) 



dtT = 7-IAT 



where 7 is the average of the flip rates ^x deflned by 

One of the main interest of the model is its non-reversibility. To the best of our knowledge, it is the first 
time that hydrodynamic limits are established for a non-reversible interacting particle system evolving in a 
random medium. In fact, our first motivation was to work with a simplified version of the energy conserving 
model of heat conduction with random masses ([2]) and we think that some of the methods developed in 
this paper could be useful to study this model. 

The derivation of the hydrodynamic limits presents three difficulties: the first is that the system is non- 
gradient. The second one is that it is non- reversible and that the symmetric part S of the generator is 
very degenerate and gives only few pieces of information on the ergodic properties of the system. The third 
difficulty is more technical. The state space is non-compact and the control of high energies is non-trivial. 
The first problem is solved by using the "corrected empirical measure" method introduced in [12], [14] and 
some special features of the model. For the second one, we apply in this context some deep ideas introduced 
in [TT] (see also [IB])- The third problem is solved by observing that the set of convex combinations of 
Gaussian measures is preserved by the dynamics. The control of large energies is then reduced to the control 
of large covariances. 

In the perspective to study heat conduction models with random masses our main interest lies in the 
properties of the diffusion coefficient (given here by I/7). 

The diffusion coefficient is also often expressed by the Green-Kubo formula, which is nothing but the 
space-time variance of the current at equilibrium. The Green-Kubo expression is only formal in the sense 
that a double limit (in space and time) has to be taken. For reversible systems, the existence is not difficult 
to establish. But for non- reversible systems even the convergence of the formula is challenging ([E]). Let us 
remark that a priori the Green-Kubo formula depends on the particular realization of the disorder. 
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If we let aside the existence problem, widely accepted heuristic arguments predict the equality between the 
diffusion coefficient defined through hydrodynamics and the diffusion coefficient defined by the Green-Kubo 
formula. 

The second main theorem of our paper shows that the homogenization effect also occurs for the Green- 
Kubo formula (see Theorem [2]): for almost every realization of the disorder, the Green-Kubo formula exists 
and is independent of the disorder. Unfortunately we did not succeed to prove that the value of the Green- 
Kubo formula is I/7. 

The paper is organized as follows. In section [2] we define the system. The proof of hydrodynamic limits 
is given in section [31 The two main technical steps which are the derivation of a one block lemma and the 
control of high energies are postponed to sections S] and [5] The study of the Green-Kubo formula is the 
content of the last section. 



2. The model 

For any a > 0, let fla be the set composed of configurations rj = {7]x)xe'i such that ||77||q < +00 where 

Let ri — C^a>Q^a be equipped with its natural product topology and its Borel cr-field. The set of Borel 
probability measures on will be denoted by 'P{^). We also introduce the set Cq (51), fc > 1, composed of 
bounded local functions on VL which are differentiable up to order k with bounded partial derivatives. 

The time evolution of the process {r]{t))t>o can be defined as follows. Let {J\fx ; a; G Z} be a sequence 
of independent Poisson processes. We shall denote by 7a; > the intensity of J\fx- We assume there exist 
positive constants 7_ and 7+ such that 

(3) Vx e Z, 7- < 7x < 7+ 

For every realization of the random element Af = {Afx)x<£Z, consider the set of integral equations: 

(4) rUt) = (-l)-^^(*) (^vx{0) - l\-l)^^^^Hvx+i{s) - Vx-M}ds^ 

For each initial condition cr e il the equations (jj]) can be solved by a classical iterative scheme. The 
solution r]{-) := rj{-, a) defines a strong Markov process with cadlag trajectories. Moreover each path rj{-, a) 
is a continuous and differentiable function of the initial data cr ( [5] , [10] , [H] ) . We define the corresponding 
semigroup (Pt)t>o by {Ptf){a) = Mj^{f{i]{t,a))) where Ej^ denotes the expectation with respect to the 
Poisson clocks and / is a bounded measurable function on fl. 

Since the state space is not compact Hille-Yosida theory can not be applied directly. Nevertheless, the 
differentiability with respect to initial conditions and stochastic calculus show that the Chapman-Kolmogorov 
equations 

{Ptf){^) = fi'y) + f\cPJ){a)ds, f G c^ift) 
Jo 

and 

(Ptm^) = fi<y) + f\PsCf){a)ds, f e c^m 

Jo 

are valid with C the formal generator defined by ([1]). 
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The two Chapman-Kolmogorov equations permit to deduce that the probabihty measures v g 'P(i^), 
which are invariant for {ri{t))t>Q, are characterized by the stationary Kohiiogorov equation 

J {CJ){j^)du{7^) = for all / e CH^l) 

In particular, every Gibbs measure with inverse temperature /3 > is a stationary probability measure. 
Observe that /i^ is nothing but the product of centered Gaussian probability measures on R with variance 
/3~^. It is easy to show that (-Pt)f>o defines a strongly continuous contraction semigroup in L^(/i^) whose 
generator is a closed extension of C. 

In fact, the infinite volume dynamics is well approximated by the finite dimensional dynamics vP{t) = 
{vi^i{t) ; X e Z}, n > 2. It is defined by the generator £„ — An + Sn where, for any function / e Cq (fi), 

n-l 
x——n-\-l 

and 

n 

a;— — n 

Observe that 77"(i), |a;| > n, do not change in time. Moreover, the total energy J^xez^x is conserved by the 
finite dimensional dynamics. We denote by (P")(>o the corresponding semigroup. Let us fix a positive time 
T > 0, a parameter a > and a function (p € Cg (fi). One can prove there exist constants C„ := C(n, a, T, 0), 
n > 2, such that 

(5) sup \{Pr4>){v)-{m{v)\<CnMl 

tG[0,T] 

and 

lim C„ = 

This approximation is only used in the proof of Lemma [9l The proof of ([5]) in a similar context can be 
found in [3], chapter 2 (see also [H]). 

3. Hydrodynamic limits 

For any function u : Z R, the discrete gradient Vu of u is the function defined on Z by 

Va; e Z, {Vu){x) = u{x + 1) - u{x) 

The hydrodynamic limits are established in a diffusive scale. This means that we perform the time 
acceleration t — N'^t and the space dilatation x — x/N. In the rest of the paper, apart from section [HI the 
process (?/(i))f>o is the Markov process defined above with this time change. The corresponding generator 
is N^C. 

The local conservation of energy = 77^/2 is expressed by the following microscopic continuity equation 

eAt) - ex{0) = -N^ f (Vj,_i,,)(77(s))ds 
Jo 

where the current jx.x+i '■= jx,x+ii'>]) is defined by 

jx.x+iiv) = -VxVx+l 

We denote by Co(IR) the space of continuous functions on K with compact support and by Cq (R), fc > 1, 
the space of compactly supported functions which are differentiable up to order k. Let M (resp. A^+) be 
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the space of Radon measures (resp. positive Radon measures) on M endowed with the weak topology. If 
G e Cg (R) and m € M then (m, G) denotes the integral of G with respect to m. 

The empirical positive Radon measure tt^ G ^ associated to the process e{t) :— {ex{t) ; x E Z}, is 
defined by 

Fix a strictly positive inverse temperature profile /?□ : IR — > (0, +cx)) and a positive constant /3 such that 

H 2 

= 



(6) lim 



JV->-oo N 



1 1 



-1. |_ 

Denote by /j,^ — ^ ^(^) the product probability measure defined by 



where gp{u)du is the centered Gaussian probability measure on M with variance (3 ^ . 
We assume that the initial state satisfies 

(7) H{^l''\^l-p) < CoN 

for a positive constant Cq independent of N. Here H{-\-) is the relative entropy, which is defined, for two 
probability measures P,Q £ Vift), by 



(8) 



H{P\Q) = sup |y HP - log (^J e^dQ^ | 



with the supremum carried over all bounded measurable functions on 51. Let us recall the entropy inequality, 
which states that for every positive constant a > and every bounded measurable function 0, 



(9) 



>dP < fl-i jlog (^J e'^^dQ^ +H{P\Q)^ 



Fix a positive time T > 0. The law of the process on the path space D{[0, T], fl), induced by the Markov 
process {v{t))t>o starting from , is denoted by P^n. For any time s > 0, the probability measure on fl 
given by the law of 7y(s) is denoted by /z^. 

Since entropy is decreasing in time, (O implies that 

(10) Vs>0, H{^^^\^if,)<CoN 

The conditions dH) and ([7]) are introduced to get some moment bounds (see section [5]) . They are satisfied 
by any continuous function going to /3~^ at infinity sufficiently fast. 

Theorem 1. Let {^x)xi£Z be a sequence of positive numbers satisfying ^ and such that 

(11) lim t;X!^^ = ^' J™ ^ H ^^ = ^ 

x—1 X——K 

for some 7 G (0, 00). Assume that the initial state = /^^(.) satisfies ^ and f3o satisfies 
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Then, under ¥^n, tt" converges in probability to Tt/2 where Tt is the unique weak solution of (0); For 
every G G Co(K), every t > Q, and every 5 > Q, 



lim I 

Af-s-oo 



(7rf,G)-i(T„G) 



> 5 







We follow the method of the "corrected empirical measure" introduced in [121 lUl ■ Since the state space 
is not compact, technical adaptations are necessary. In particular, it is not given for free that the corrected 
empirical measure and the empirical measure have the same limit points for the weak convergence. It would 
be trivial if the state space was compact. Moreover, a replacement lemma, reduced to a one-block estimate, 
has to be established (see section U]). 

For any G € Co(IR), we define T^G : Z ^ M by 



{T,G){x)=Y.^rj+lj+i)\G 



j<x 




Observe that 



N- 



1 



■ [{T^G){x + 1) - {T^G){x)] = {\/nG){x/N) 



Ix + lx+1 

where stands for the discrete derivative: {VnG){x/N) = N{G{ix + l)/N) - G{x/N)}. 

Since T-yG may not belong to £i(Z), we modify T^G in order to integrate it with respect to the empirical 
measure. Fix < 6* < 1/2 and consider a increasing nonnegative function g defined on M such that 
g{q) = for g < 0, g{q) = 1 for g > 1 and g{q) = q foT q G [9,1 - 9]. 

Fix an arbitrary integer i > and let 5 = ge,^ : K — >■ R be given by 



giq) = g{q/£) 



We define 



{T^,iG){x) = {T^G){x) 



7,G 



iTjg){x) 



where 



T^,h - J2(lx + Ix+i) {h{{x + l)/N) ~ h{x/N)} 



In the rest of the paper we make the choice t :— i{N) — N^^'^. 



Lemma 1. For each function G E Cq( 



id each environment 7 satisfying and ill}) . 



lim A^i/" sup \T^,iG{x) - 7G(a;/A^)| = 



N—^oa 



xel 



and 



lim N^/^T^,G = 



Ar->oo 



Proof. This is a slight modification of Lemma 4.1 in [14] . 

We shall denote by e A4 the corrected empirical measure defined by 

(G) = X^'-^iG) = l^r^,,G(x)e,(x) 



□ 



xei 



7 



The system is non-gradient but we have 



Jx,x+1 — 



7x + 7x+i 



ex + ^x-iVx+i 



1 



2(7x +7x+i) 



VxVx+l 



This imphes that 



{AnG){x/N) - ^{A^g){a:/N) 



ex + -jTjx-irix+l 



(VivG)(a;/7V) - ^{yNg){x/N) 



VxTlx+l 



where Ajv stands for the discrete Laplacian: 

{AnG){x/N) = {G{{x + 1)/N) + G{{x - 1)/N) - 2G{x/N)} 

Therefore, we have 

(12) (G) - Xo^(G) = C/f (G) + Vf{G) + (G) 

with M^{G) a martingale and t/^(G), V^{G), which are given by 



where 



and 



i7f (G) = fds ^Y,B^{x/N) ( e,(s) + ^,.^{s)vx+i{s) 

x& \ ' 

B%{xlf^) 



{AnG){xIN) - '^{ANg){x/N) 



{VnGKx/N) - ^(VnoKx/N) 



'■1,9 



{r]j,{t)f]x+l{t) - r]a:{0)f]x+l{0)) 



Lemma 2. The sequence {{X^J^ n^ds) e D{[0,T],M) x D{[0,T],M+) ; N>l} is tight. 
Proof. It is well known that the sequence 

Trfds) e Di[0,T],M) X D{[0,T],M+); N > 1^ 



is tight if and only if the sequence 



X (G), / Tr^{H)ds] e D{[0,T],R)x D{[0,T],R); N>1 



is tight for every G,H € G^(M). 

By Aldous criterion for tightness in D([0, T], M)^, it is sufficient to show that 
(1) For every t G [0, T] and every £ > 0, there exists a finite constant ^ > such that 

supP^^(|Yj^(G)| > A) <£ 



(2) For every S > 0, 

where 6 is the set of all stopping times bounded by T. 



lim lim sup sup F^n \ I Y^^.g {G) - {G)\ > 5] = 

e^O Tv^oo rGe,6><e 
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for r^(G) = X^{G) and Y.^ (G) = /„ (G)ds. 

Since G has compact support, there exists a constant K > (independent of t and N) such that 



(13) 

and consequently 
(14) 



< ( N^/Uvip\jG{x./N) - {T^^iG){x)\) V 



« [|7(7rf,G)-Xf(G)| 
1 



Ar5/4 



\x\<KN^/'^ 



rfi 7(<,G)~Xf(G) 



< (iVi/4sup|7G(x./iV)-(r^,,G)(x)| 



dtE„ 



E 



a;|</^JV5/4 



By Lemma [T] and Lemma [TUl the right-hand side of (resp. of (HU) vanishes as — > cxi. Hence, it is 
sufficient to show Aldous criterion for Y.^ {G) = X^{G) and for r^(G) = (G)ds. 

From the definition of the Skorohod topology, it is easy to show that the application $ from D([0,r],M) 
onto itself defined by 

$ : X := {x{t) ;0<t<T}^ <^{x) |^ x{s)ds ; < t < r| 

is continuous. Thus, if {X^{G))n is tight, then {J^ {G))n is tight. 
Therefore it just remains to show Aldous criterion for Y^{G) = X^ (G). 
Proof of (1) for X.^(G); 



l^f (G)| >A]< jE^. ( 1e \iT,,iG){x)\ e.(i) ) 



We write {T^^eG){x) = {{T^jG){x) - ^G{x/N)) + jG{x/N) and we get that 



[|Xf(G)|>A] < ^IEp«|^ E l(T7/G)(x)-7G(x/A^)|e,(i) 



N 



\x\<KN^/* 



I ^ E \G{x/N)\eS) 



N 



\x\<KN 



The first term on the right-hand side of the previous inequality can be bounded above by the right-hand 
side of (IT51) . which vanishes. By Lemma [TUl the second term is bounded above by G/A with a constant G 
independent of N. Therefore, the first condition is satisfied. 
Proof of (2) for X.^(G): 

Recall the decomposition p^. In order to estimate the term 



[\U^^^{G)-UAG)\] 
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we observe that \B'^{x / N)\ is bounded above by 



C 



^\x\<KN 



\2\g\ 

IT.. J 



(£-2 + (7Vr3)) 



6[i-2e4+26»]u[-2e,2e]} 



where C, K are constants depending on 6 and G but not on N . By Schwarz inequahty, we are reduced to 
estimate 



and 



IT,,, I \P Nl\ 



(is 



N 



|2;|<(l+3e)Af£ 



By Lemma [TUl the first term is of order one. It is not difficuit to show that iiminfjv-!.oc T~f_g > 0, and Lemma 
[T] gives iTj^G 0. Thus, by Lemma [TOl the second one vanishes as N goes to infinity. 
The two last terms of ([T2|) are given by 



(VjvG)(x/iV) 



Ty,G 



= y,^(G)+Aff(G) 
By using Lemma [1] and Lemma I10[ similar estimates as before show that 

lim sup E^iv [\V/'{G)\]^0 

te[0,T+e] 



VxVx + 1 (s) 



By computing the quadratic variation of the martingale (G), one obtains that (we recall that sup^ jx < 



7+) 



(Aff+,(G)-Aff(G))^ 



< 



274 
N 



2 



1:G 



{Vn9){x/N) ) es{x)es{x + l) 



7:9 



Observe that 



{VnG){xIN) 



7:G 



7,9 



Since lim inf jv-s-oo Tj^g > and iT^.c ~^ 0, from Lemma [TOl we get 



s>0 



\x\<KN 



and 



SUpE^« l—^ ^ 6^(4 

go to with N (and are in particular bounded above by a constant independent of N). 



□ 
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Lemma 3. Let (a,/3) M x Ai^ be a limit point of the sequence 

X.^,^ Trfds^ e D{[0,T],M) X D{[0,T], M + ) ; N > 1 
For every G G C'o(M) and every t € [0,T], we have 
at(G)-ao(G) =7-1 
Proof. In the proof of the tightness of we have seen that the term 
E 



,(AG)ds, Pt^r 



VxVx+i (s) 



and the term 



\ T Q ( 1 \ 



vanish as iV ^ oo. By using Lemma [5l it implies that 

at(G)-ao(G) = /3t(AG) 

Moreover, by (|14p . we have 



□ 

Lemma 4. ^nj/ /imit point /? 0/ f/ie sequence {J^ir^ds G £'([0, T], A^+) ; > 1} is suc/i t/iai, /or any 
t G [0,T], /3t is absolutely continuous with respect to the Lebesgue measure on M. 

Proof. Fix a positive time t and let R^n be the probabihty measure on given by 



for every Borel subset A of + . Let J : [0, +00) be a continuous and bounded function. By the 

entropy inequaUty ([9]) and by using (jlOp we have 



(15) 



J(^)di?^«(7r) < Go + ^log e'^'^-UR^^i^)^ 



By the Laplace- Varadhan theorem, the second term on the right hand side converges as N goes to infinity 
to 

sup [J(7r) - Jo(7r)] 

where /q is the large deviations rate function for the random measure tt under R^^-^ . It is a simple exercise 
to compute the rate function Jq. We have 

/o(7r) = sup 

fee am 

where Mp(a) is the Laplace transform of ?7q/2 under /i^: 



f{u)n{du)~ / logM^-(/(u))d7. 
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if a < /3, and +00 otherwise. 

The function /o also takes the simple form 



h{'!T{u))du if Tr(du) — TT{u)du, 
+00 otherwise 



where the Legendre transform h of Mp is given by h{a) = f3a — 1/2 — l/21og(2a/3) > if a > 0, and +00 
otherwise. 

Let {fk)k>i be a dense sequence in Co(M) with /i being the function identically equal to 0. Then Iq is 
the increasing limit of Jk > defined by 



sup 

l<J<k 



fj{u)^{du) - / log M-p{fj{u))du \ A k 



By using (fTS]) we have 



limsup / Jfc(7r)(ii?^jv (tt) < Cq 
for each k. Since Jk is a lower semi-continuous function, any limit point R* of is such that 

Jfe(7r)di?*(7r) < Co 



By the monotone convergence theorem, we have / Io{TT)dR* (tt) < Cq < +00. Since Io{tt) is equal to +00 if 
TT is not absolutely continuous with respect to the Lebesgue measure, it implies that 

R* {7r;7r((iu) = Tr(u)du} = 1 

and the lemma is proved. □ 

We conclude as follows. Let (a, (3) be a limit point of (X.^(G), Jg (G))n>i- From the equation 



a.(G)-ao(G) = 7- 



asiAG)ds 



we see that a is time continuous. Moreover, if ^ is a subset of R with zero Lebesgue measure, then /3t{A) — 
for any t G [0,T]. This implies that at{A) = for any t £ [0,r], i.e. that at is absolutely continuous with 
respect to the Lebesgue measure on R. 

By uniqueness of weak solution to the heat equation, we have that 2a is the Dirac mass concentrated 
on the (smooth) solution of the heat equation {t,u) £ [0,T] x R — > Ti(it) starting from j(3q^: dtTt = 
7-iAfi, fo = 7/3o-'. 

Hence we conclude that {X!^ £ D{[0,T], A4) ; > 1} converges in distribution to (T.(u)/2) du. Since the 
limit is continuous in time we have that {X^ X > 1} converges in distribution to the deterministic limit 
(Tf (m)/2) du. Since convergence in distribution to a deterministic variable implies convergence in probability, 
this implies that 

1 



lim 

N^oc 



X^{G) 



Tt{u)G{u)du 



We use again and the fact that 7 ^Tt = Tt to get 



lim I 

N^oc 



<{G)--l Tt{u)G{u)du 



> £ 



> £ 







= 



and the theorem is proved. 
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4. One-block estimate 
The aim of this section is to prove the following so-called one block estimate ([15j). 
Lemma 5 (One block estimate). For any G £ Cq{M.), any t > 0, and any S > 0, 



lim 



1 /■* 

-^(AjvG)(x/7V) / 77.-1 (s)?7,+i(s) 



> S 



= 



Since G G C^(M), we can replace {AnG){x/N) by 



2k 



- {A^GXy/N) 



\y — x\<.k 



as soon as fc <C and we are left to prove that 

1 



lim lim E„ 

fe— >oo N—^oo 



2N- 



|.|<A 



ds 



2k 



\x — y\<k 



Given two probability measures P, Q on 17 and A a finite subset of Z, Ha{P\Q) denotes the relative entropy 
of the projection of P on with respect to the projection of Q on R'^. We shall denote the projection of 
P on by P|a- If A = A^ = {— fc, . . . , k}, we use the short notation P^. 



We define the space-time average of {Hs)o<s<t by 



1 



{2N + l)t 



E 



\<N 



TxfJ.^ ds 



Here denotes the shift by x: for any 77 G 51, the configuration T2:77 is defined by {Txffjz = Vx+z', for 
any function g on f2, T^g is the function on given by {Txg){ri) = g{Txri); for any p G 'P{^), r^p is the 
push- forward of p by Ta; . The probability measure p is said to be translation invariant if r^p — p for any 

xeZ. 

We have to show 



(16) 



lim lim 



) A— >C« ImA 



dv, 



N 



1 



2k 



1 E ^y^y+'^ 



\y\<k 



= 



Lemma 6. For each fixed k, the sequence of probability measure {v^)N>k on MJ^'' is tight. 
Proof. It is enough to prove that there exists a constant Cfe < 00 independent of N such that 



(17) 

We begin to prove that 
(18) 



HA,i:^''\l^f,) = H[,y^\f,^\A,) <Co\Ak\ 



Fix a bounded measurable function (p depending only on the sites in A := A/j = {— fc, . . . , k}. Assume 
for simplicity that 27V -f 1 = {2k + l){2p + 1) for some p > 1. Then we can index the elements of the set 
{—N, . . . , N} in the following way 

{-N, . ..,N} ^ {xj +y;j = ~p,. ..,p; y e A^} 
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where Xj = 2kj + 1. Since (f> depends only on the sites in A, it is clear that under /i^, for each j/ G A, the 



random variables [Tx^+yip) ^ are i.i.d. 



Let fjL^ = t ^ /g ds. By convexity of the entropy and p^ . we have H{fif\fi^) < CqN. 
We write 



< 




2N 

< 



where we used the entropy inequality ([9]) and the convexity of the application / — > log (/ d/i^e-^). By 
independence, for each y, of {Txj+y<f) ^ and the translation invariance of /i^, we get 

(19) j cj,dv^ < Co\A\ + log (^J e'^'d^l^^ 

This implies (|18p and, by the entropy inequality ([9]), the inequality ([TTl) . □ 

For any k, let be a limit point of the sequence {v^)N>k- The sequence of probability measures {v^)k>o 
forms a consistent family and, by Kolmogorov theorem, there exists a unique probability measure v on 
such that Vk — i^k- construction, the probability measure v is invariant by translations. 

Lemma 7. There exists Cq such that for any box Ak — {—k, . . . , k}, k > 0, 

(20) HAMf^p) < Co\Ak\ 
Proof. We have seen in the proof of the previous lemma that 

(21) i?A.(^'^lM/3) = H (^.f I /i^U,) < Co\Ak\ 
Since the entropy is lower semicontinuous, it follows that 

□ 

A translation invariant probability measure z/ on f2 such that (j20l) is satisfied is said to have a finite 
entropy density. By a super-additivity argument (see [3], [H]), the following limit 

(22) H{v\^ip) = hm — 

fc^oo |Afe| 

exists and is finite. For any bounded local measurable function </) on fJ, we define the limit 



14 



CEDRIC BERNARDIN 



The entropy density H{v\ijLp) can be expressed by the variational formula 
(23) H{iy\n0) = sup j [ cpdiy - F{(j3) 



where the supremum is taken over all bounded local measurable functions on fi. 
We now show the following lemma 



Lemma 8. For any function F £ Cq{VI), we have 



CFdiy = 



Proof. Assume that F £ C^{^}) has a support included in R'^'=-i. We have 



CFdi^ 



£F dvk — lini / CF dvi, 



N 



Define G = {2N + E|x|<iv ^^F- By Ito formula 



N- 



d^,^ ir,)G{r,) - d^i^ir,)Gi7^) 



= J^dsj di,^{r^){CG){r^) 
j dfl^iv) {CG)iv) 
J di^^irj) iCF)irj) 

Since F (and hence G) is bounded, the left-hand side goes to as iV goes to infinity and it follows that 

CFdiy = 



□ 



Recall that we want to show (1161) . From the previous lemmas, it is sufficient to prove that 



\y\<k 



VyVy+l 



for any p £ Vi^) such that p has finite entropy density, is stationary for C and translation invariant. 

Proposition [T] gives the characterization of stationary probability measures, translation invariant, and 
with finite entropy density. By using the notations of this proposition, to complete the proof of Lemma [Sj 
we have to show that 



lim / dA(/3) 

'^-^ac 7(0, + oo) 



1 



2k- 



\y\<k 
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lim / d/ii 



Since, under /i^, the random variables {\/^r]y)y are distributed according to standard independent Gaussian 
variables, and J I3~^d\{l3) < +cxd, it remains to prove 

\v\<k 

By using Scliwarz inequality, a simple computation gives the result. 

Proposition 1. Let v he an invariant measure for C which is translation invariant with finite entropy 
density. Then, v is a mixture of the Gaussian product measures fip, /3 > 0, 

(0,+oo) 

and the probability measure A on (0, +oo) is such that 

l3~^dX{l3) < +0O 



'(0,+oo) 

In order to give the proof of this proposition, we need the following lemma 

Lemma 9. Let v he an invariant measure for C, translation invariant with finite entropy density. Then, 
for any local measurable bounded function (p on fl, we have 

Vx e Z, J [<j){r]^) - (t>{7j)] di^^O 

Proof. We only give a sketch of the proof since the arguments are almost the same as in [TT], Proposition 
6.1 (see also chapter 2 of [3]). 

The proof is divided in two steps. Let us first consider a generic probability measure i^,, not necessarily 
translation invariant, such that < +oo and let us denote by g the density of i^* with respect to /^^. 

We introduce, for any n, the Dirichlet forms 

(24) Dn{iy,)^ sup I - y ^du, | 

where the supremum is carried over the set T composed of the positive functions : — (0, +oo) such that 
< <il)<MioY some positive constant M . 
It is easy to check that if Dn{v*) < +oo then 



(25) 



1 " r 



where for any function u : — > M, YxU is the function defined by {Yxu){ri) — u{ri^) — u(ri). Observe that 
= ~2Y, so that -Sn = (1/2) E"=-„ 7x>;'- 

In fact, even if i/.^ is not absolutely continuous with respect to /i^, the Dirichlet form Dn{v*) defined by 
([24| makes sense in [0, +oo]. 

Recall that {P^ )t>Q is the semigroup generated by the finite dimensional dynamics introduced in section 
[2l We have the following well known entropy production bound (see [3j or [15], Theorem 9.2) 

Hiv^Pl' I + tZ?„(<t) < H{v, I ^i-p) 

where = i^^P^ds. 
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Let US denote the density of z/"( with respect to fi^ by g". Since < +00, we have Dn{v^^t) < +00 

and, by the exphcit formula (|25p of the Dirichlet form, 



(26) H{v.P^ I + 2J ^ j {y, ' d^lp < H{v. I 

The second term on the left-hand side of the previous inequality is composed by a sum of positive parts. We 
can restrict this for any m < n. By using ([SJ and the variational formula for the Dirichlet form, we get 
that, for any function ip G CQ{n) and any functions ipj G J^, G {—rn, . . . , m}. 



j — — m 



Pl^<j,du. - log j eUf,^ + V ^ - 
We let n — > 00 and, by JS]), we have 
(27) 



. We recall that denotes 



j Pt(j)dv^ - log / e*'d//^ + 2^ -T — diy^.t < H{v^\^i^) 

j — — m 

where P^^^t = Jq v*Psds. 

In the second step of the proof we apply ([27l) to z^* = v'^^ = v ® l^'p 

the box {— m, . . . , m} and Af„ stands for Z\A„j. Observe that H{v^^\^ip) = H\^^^{iy\^p) so that 

lim {2m+l)-^H{iyi"''>\ng) = H{iy\fig) 

By choosing (p = X^Il-m '^^'/'o, V'l = '^iV'o, with (po e C(^(ri) and V'o e -7^, we get 

^O^V'O ,_(m) ^ rr/ (m)| 



t—~m 

We claim that 



(28) 



lim — - [ Pt{n(l)o)diyi"'^ ^ ( Pic^^dv ^ ( 4>odu 

woo 2m + 1 . ^ J J J 



hm j Ti— — dvl ^' = / — ; — dv 



2m + 1 . ^ J ■00 ' 7 "00 

I— — m 

Then, by using ([23]) and optimizing over (po and V'o, we get 

f ^0'^°^ n 
sup / — ; — ai^ = u 

It is clear that we can repeat the argument substituting Yj to Yq, and we obtain 

sup / — — di' = 



so that, by summing over j, we have Dn{v) = which implies that v is invariant by any flip. 

It remains to show (pS)) . The difficulty comes from the fact that even if the function u is local Ptu is not. 
But it is easy to see, by using ([5]), that we can replace the semigroup of the infinite dynamics Pt by the 
semigroup of the finite dimensional dynamics P", if n is sufficiently large. The function P"u is then local 
and the ergodic theorem permits to conclude. 
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We refer the interested reader to [3] for the details of the arguments. 

□ 

Proof of Proposition]^ By Lemma [SJ we have / Sg dv — for any bounded measurable function g on 57. It 
follows that for any g £ C^{fl), 



(29) J Agdv^O 

Since v has finite entropy density, we have J eodv < +oo. By translation invariance, the ergodic theorem 
gives the existence i/ a.s., and in L^(z^), of 



u(ri) = lim — T]x, Eiv) = lim — rii 



\x\<e ^' ' ' \x\<i 



Since v is invariant with respect to any flip, we have z/ almost surely that uiji) — 0. 

Assume first that is exchangeable. 

For any z G [0, oo) let Vz. be the probability measure 

If z = then is the Dirac mass concentrated on the configuration Sq with each coordinate equal to 0. 

Let us now assume that z 7^ 0. 

Consider a test function g in ([29| of the form 



giv) = fiv) X 



with /, X compactly supported and smooth. It is easy to show, by taking the limit ^ — !■ 00 in (|29p with g as 
above, that 



Af dvz = 

and this can be extended to any / g Cg(r2). We apply the previous equality with a function / of the form 

f(ri) ^ VxHv) 
with (j) G CQ(n) independent of ri^. Then we get 

= JiVx+i ~ Vx~i)4>di^z{il) + X! y(%+i ~ Vv-i)Vxdjj^(j)di^2^ 

y^x 

iVx+i - ■r]x^i)(l)dvz.{r]) 



+ I iVx - rix~2)rixdn^^-,4)diy:^ + J {rix+2 - Tlx)r]xdr,^+i4>d'^2i 

J ^'^y+l ~ Vv-l)Vxdrj^(j)dVy. 



y^x — l,x,x-\-l ' 



We claim that the last term is equal to zero. This is a consequence of the exchangeability of v^.- Let A be 
the support of (f> (which does not contain x by assumption) . Observe that , for any y ^ x— l,a;,a; + l, the site x 
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does not belong to the support of {l]y-^-l — l]y^l)^■,J^(f>. Let t be sufficiently large (e.g. t > \x\+maxs^A |s| + 10). 
By exchangeability we have, for any fc > 0, that 



Hence, we get 



J iVy+i - Vy-i)Vxdr,^(f>diy^ = J {ijy+i - 'ny-i)vt+kd^^(t)diy.. 



k=0 ' 

Let £ go to infinity and use the convergence of X)fe=o to u{ri) = to conclude. 
The same argument shows that 

j 'n^^2'nxdr^^_^4>dv^ =0, y ?7x+2?7x5^,+i0rf;^z 

and, similarly, we have 

vldr^^-i4>dvz =z dr^^_^(j)dv^, I rildr^^^^(j)dv^ = z I dri^^^cjidv^ 



Hence, we proved that, for any x e Z and for any function (/> e Cq{VI) such that x does not belong to the 
support of (/), 

{rix+i - ■rix-i)4'dv^{ri) + z J {d,j^_-^ - a^^+J^dz^z = 

We apply this for a function (j) depending only on {'q2k)kei, so that, for any fc, (p is independent of i]2k+i- 
We have 



/ 



iV2k+2 - V2k)(t>dv2.{v) + Z / i9r,2k - dr,2k+2)(t>dl^z = 



This implies that the law of {rj2k)k<=:Z under is a product of centered Gaussian probability measures on M 
with variance z (see e.g. [13]). 

The same result occurs for the law of (772fc+i)fcez. 

Let now ^Aiv) — nseA'^s(^«) ^ ^^^^ function with A a finite arbitrary set of Z and real valued 
bounded functions. We write the set A in the form U Ai where is the set composed of elements of 
A which are even and Ai the set composed of elements of A which are odd. Let i?o be a set composed of 
even sites such that |i?o| = |^i| and AC\ Bq — %. Let ct be a permutation on Z such that (j{Ai) = Bq and 
Aq is fixed under the action of u. We denote by ct • ry the configuration defined by (ct • ri)x = rj^^x)- By 
exchangeability of we have 

and 

'i>A{cr-v)= n n 'l^<^-Hs){Vs) 

Since the function $yi(f - J?) is a function depending only on {f]2k)k£i and Ao H i?o = we know that 
v^{^a{(J-v)) = n ( / 't^s{x)gi/^{x)dx\ ]J f / (l)a-Hs){x)9i/z{x)dx 

We recall that g^ji, is the density of the centered Gaussian probability measure on R with variance z. Hence, 
we proved 
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which shows that dv2,{rf) is equal to Y[x£z9i/2iVx)dr]x- 

We now show that ly is exchangeable. Let us consider the test function xiv) — 't'iVx, Vx+i)''P{^yi V 7^ -'c, x+1) 
with (j), tp smooth and compactly supported functions. By (j29p . we have 



duAx^O = J A(j)ijdi^ + J (j)Aipdv 
Observe that the second term is given by 

X : Vx+1 ) 

y^x,x+l 

This is equal to zero because v is invariant by the flips and the function r] — > r]y{de ^j){vi){riy^i—'qy^i)(j){rjx, r)x+i) 
is an odd function of ?7j^ iov y ^ x,x + 1. 
Moreover we have that 

{M){->l) = {Vx+2 - ■nx)dr,^^,4> + iVx+l - ?7x-l)<9,,> 

Remark that ?ya;+2'0f^77a;+i0 is odd with respect to rix+2 so that its integral with respect to v is equal to 0, 
and similarly for rix-iipdrf^(t>- Hence, we get 

dvi-q) {r]a;+id,j^(j) - rixd„^^^(f>) ip = 

This equation implies that i>{rix,rix+i\{ey]y ^ x,x + 1)) is exchangeable. 
Let now $ be a local test function of the form 

where {4>s)s is a sequence of bounded smooth functions equal to 1 for \s\ > A for a positive constant A. Our 
aim is to prove that for any x we have 

(30) = K<J>(^)) 

which implies the exchangeability of v. We can assume that each 0^ is even or odd since every function can 
be decomposed as the sum of an even and an odd function. Moreover each even function (j)s(r]s) takes the 
form (psiss) for a suitable function (jj^. 

If one of the (f>s is odd, since v is invariant by all flip operators, ([30l) is trivial because the two terms are 
equal to zero. We assume that all the (ps are even so that $ is in fact a function depending only of the 
energies Cg and we write = $(e) = Osez 'Psies)- We shall denote by v the law of e := {cy ; y G Z}. We 
have 



^{ri)diy{ri) = J di>{e)^{e) 

di){ey;y x,x + l) ^{e)di>{ex, e^+ilcy, y ^x,x + l) 
dv{ey]y x,x + l) (^J ^{e'^'''^^)di'{ex,ex+i\ey,y ^ x,x + 1) 

x,x+l\ 



where we used the exchangeability of (^(r/^;, ?72;+i|(ey; ?/ ^ x,x + 1)) in the third equality. It concludes the 
proof that v is exchangeable. 
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Hence, we can express i' as a, mixture of /i^, /3 G (0, +00], with the convention that /^oo is the Dirac mass 
concentrated on the configuration So- 



'(0,+oo] 

In fact, A is the law under u of the random variable l/£{'q). It remains to prove that i'{£{ri) = 0) = 
A({+cx)}) = 0. It is a simple consequence of the fact that Ha^.^i'Iii^) < Co|Afc| for any k and in particular 
for fc = 0. By ([8]), we have that for any positive real M 

Co>M J l[o}{x)d,y\^^^{x)~log (^J e^'^i'>y^^'>g-^{x)dx^ = MA({+^}) 

Since M is arbitrary large, it follows that A({+oo}) = 0. 



□ 



5. Moments bounds 



The aim of this section is to give the proof of the following lemma: 

Lemma 10. Let jj,^ be the probability measure associated to a temperature profile bounded below by 

a strictly positive constant such that (0j and ^ are valid. Let (Mjv)7v>i be 0, sequence of positive integers 



such that liminfjv- 
that 



, Mjq /N > 0. Then, there exists a positive constant C , which is independent of N , such 



sup 

t>o 



1 

Mj7 



E 

|2:|<A/jv 



ex{t) 



< C 



lim sup 



>oo t>0 



Ml, 



E 



elit) 



Let us first explain why the second equality of this lemma is nontrivial. The standard arguments to get 
moment upper bounds are based on the entropy inequality ^ and the existence of exponential moments. 
In our case it would be necessary to have /i^(e"''o) < +00 for a sufficiently small. This is false since fip is a 
Gaussian measure. In [T], following an idea of Varadhan, and despite the absence of exponential moments, 
the use of the entropy inequality for the microcanonical measure was sufficient to get a weak form of the 
lemma we want to prove. This approach cannot be carried here because we are in infinite volume and because 
the Dirichlet form is too degenerate to reproduce the argument. 

Proof. The first statement is a simple consequence of the entropy inequality Indeed, for any S > 0, we 
have 



E„ 



1 

Mn 



E 

\x\<Mn 



< 



1 



SM 



N 



SM 



-loa 



N 



-Vl/2 



The first term on the right-hand side is of order one by ([T0|) and the second term is also of order one if S is 
sufficiently small. Hence the left-hand side is of order one in N uniformly in time. 

The bound on the second moment of the energy is more difficult to obtain and the entropy inequality is 
not sufficient. We exploit here the Gaussian structure of the initial state. 

Recall the integral equations (|4]) defining the dynamics. Each Poisson process Mx is interpreted as a clock 
and a jump of Afx as a ring of the clock. Conditionally to the realization oi Af = {Afx)x, the dynamics 
is linear, thus the law remains Gaussian in the time interval between two successive rings. When a clock 
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rings the flip operation conserves the Gaussian property of the state. Hence, conditionally to N , the state 
remains Gaussian for any time. It follows that the law fi^ of the process at time i is a convex combination 
of Gaussian measures Gm,c with mean m £ and correlation matrix C S iSz(R), the space of symmetric 
matrices indexed by Z: 

/^f = / dpt{m,C)Grn,c 



Moreover, the convex combination pt{m,C) is the law at time t of the Markov process {■m{t),C{t)) with 
formal generator N'^Q where 

{QF){m, C) — ^(Ca;+i,j, — Cx-i,a + Cx,y+i — Cx,y)dc^_yF 
x,y 

+ - mx)d^^F + Y,[F{C'',mn - F{C, m)] 

X X 

with given by 

^^a;^ _ I Cu^v if [u 7^ a; and v ^ x] ot [u = v ~ x], 
1 —Cu,v otherwise 

and 

(to-)„ = (-1)^«(--")to„ 
In other words, (C(-), are the solutions of the following integral equations 

- [Cx+i^s) - + Cx^y+,{s) - Cx,y-i{s)] ds) 

mx{t') - (m,(0)-/j'(-l)^^(*') [m,+i(s)-m,_i(s)]ds) 

with initial conditions 

m,(0) = 0, a,y(0) = Jo(a; - y)P^\x/N) 

and = tN"^. 

The existence and uniqueness of solutions is easily established (by the same methods as presented in 
section [2]) in the space H = x Hi, where 



Ho - Pi J m e ; ^ e-^l^'m^ < +oo I 

Q>0 I X J 

Hi = fl I C e 5z(R) ; e^^^l^l+I^I^C^ 

a>0 I x,y 



,<+oo| 

Observe that the initial condition belongs to H. Moreover, for any (m, C) G H, the Gaussian measure with 
mean m and correlation matrix C is meaningful (see e.g. chapter 2 of [5]). 
This Markov process conserves the three quantities 

(31) E"^- E E^--- 

The initial condition /i^ is such that po is the Dirac mass concentrated on 



= 0, C,,j, = <5o(a:-y)/3o"'(a^/iV) 
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Therefore, we have m{t) = for any t > 0. By denoting, with abuse of notations, by pt{C) the law of C{t) 
at time we have by the two last conservation laws pip that 



dpt{C) 



1 



1 



Moreover, we have 



IE,, 



|x|<M„ 



= NI 



N 



Mn 



J2 J dptiQGo^cirit) 
J2 Jdpt{C)Cl, 



\x\<Mm 

3 



|a;|<A/„ ■ 



" \x\<Mk ' N \x\<Mk ) \ / 

where we used the fact that, for a Gaussian centered variable, the fourth moment is given by three times 
the square of the second one. 
Observe that 



1 



^ \x\<Mn 



and this term is order by the first part of the lemma. 

Up to terms of order M^^, we are left with 



dptiC) 



M 



N 



E iCx..x-r')' 



x|<Afiv 



< J dptiC) 
= [ dpoiC) 



^ E {Cx.^y-r'Soix-y)f 



1 



E iCx.,y~r'So{x-y)f 



M 

since the penultimate sum is conserved by {C{t))t>o. By the assumption ([6]), the last term goes to zero as 
N goes to infinity. □ 

6. Green-Kubo formula 

In this section we study the homogenization properties for the diffusion coefficient in the linear response 
theory framework. To present the results we have to introduce some notations. 

Let {'^x)xi£Z be a sequence of i.i.d. positive random variables satisfying the assumption 

P [7- < 7x < 7+] = 1 
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where P is the probabihty measure on given by the law of the disorder 7 — {'^x)x£Z- The corresponding 
expectation is denoted by E. 

In this section, time is not accelerated by a factor N'^ . We first consider the closed system of length > 1 
with periodic boundary conditions. Let = {0, . . . , — 1} be the usual discrete torus of length N. The 
generator of the system is given by dU) with the sums over x & Z replaced hy x &Tn- 

Linear response theory predicts that the diffusion coefficient D := D{{^}, /3) appearing in ^ is given by 

(32) D = lim lim La, (A) 

where Ln '■= Lj^^ is the Laplace transform of the current-current correlation function. It is defined for 
z e H+, H+ = {zeC; m{z) > 0}, by 

?2 



x£Tn l/6Tjv 

Here, (•, •) := (•, •)^ denotes the scalar product in L^(^^) where 

is the Gibbs equilibrium measure with inverse temperature /3 > on R^". We also use the short notation 
:— (•) for the expectation with respect to /i^. 
The Laplace transform Ln can be written as 



Observe that the definition p2p is only formal since it is not clear a priori that the limits exist. 

We also consider the homogenized Green-Kubo formula for the infinite volume dynamics. It is defined by 

(33) ^(/3)= hin L^(A) 

A>0,A-i.O 

where L := is the Laplace transform of the averaged current-current correlation function. It is defined 
for z e by 

L{z) ^YJ„ ^ Joa(i),jo,i(0) » 

where <C •, • ^=<C •, • is the inner product defined for bounded local functions / and g by 



«/,ff»0=E(^[(T,/,g) 



P - {f)p{9)fi] 



We shall denote by L^(^ • ^) the Hilbert space generated by the set of bounded local functions and the 
inner product ^ •, • 

The aim of this section is to show the following homogenization result 

Theorem 2. For almost every realization 0} the disorder ^ , the Green-Kubo formulas \3l3jl and \SS\) converge 
and are equal: £'({7},/3) = D{j3). Moreover, D is independent oj 13. 

We recall that the functions Lm and L are analytical functions on (see e.g. [21], Theorem VIII. 2). 
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Lemma 11. There exists a constant C :— C(/3,7+), independent of N , 7 and z € such that 

\Ln{z)\ < C 

Proof. The proof is a simple consequence of Proposition 6.1 in [15] and of the fact that Sjx,x+i — ^'^{"fx + 
lx+i)jx,x+i (see also Theorem 2 in [2]). □ 

Let {rj; /3, 7) be the solution of the resolvent equation in L^((-)): 

{z - £N)h^ = ^ jx,x+i 

We have 

/32 / 1 \ 
Ln{z) = yK'^z 2^ jv,v+i I 

Let hz '■— hz{ri; (3) be the solution of the resolvent equation in L^(<C • 

(z - C)hz = jo,i 

We have 

L{z) = ^ < hz,jo,i > 

Observe that if ?/ is distributed according to iip then P^/'^rj is distributed according to ^i. Since hz{ri] 1) = 
hz{r}', /3) and jx,x+i is an homogeneous function of degree two in 77, it follows that L^{z) = L^(z). This implies 
the independence of the diffusion coefficient with respect to /?. 

In the following lemma we give an explicit formula for L{z) if ^^{z) is sufficiently large. 

We shall denote by P_r.vf. the law of the two-dimensional simple symmetric random walk (Sj)j>o = 
{Sj,Sj )j>o starting from (0, 1) and by Efl.vy. the corresponding expectation. Let E be the annealed expec- 
tation EEfl.vi/. . 

For any path {S'j}{j=o k} of length k, we define e{{S}k) = Y\'jZo{{Sj+i — Sj) ■ w) e {±1}, where w 

is the vector (1,1) and x • y denotes the usual scalar product of the two vectors x and y of M^. We also 
introduce the random potential 

exp{-Vz{x,y)) ^ —— ■ 

Z + Ix^yhx + ly) 

Lemma 12. There exists Aq > such that, for any z G i7+ with £H(z) > Aq, the Laplace transform L{z) is 
given by 

(34) L{z) = --Y^i-A^E [s{{S},)e-^W^^'^^SMSl - Sl) 

k=Q 

Proof. Since the generator C maps a polynomial function to a polynomial function of the same degree, the 
solution of the resolvent equation is expected to be of the form 

where 4>z{x,y), (x,y) G Z^, is the (symmetric) solution of 

5i{x)5o{y) + 5o{x)5i{y) 



(35) {z + {nx+ly)'^x^y)(iiz{x,y) + {V(j)z){x,y) = 
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with, for any function u : — > R, 

(Vm)(x, y) = {u{x, y + 1) - u{x, y-l)) + {u{x + 1, y) - u{x - 1, y)) 

We shall denote by A the real part of z € H^. In the sequel we show that, if A is sufficiently large, a 
solution to ([55)) exists, so that is of the form given above. In fact, it is not difficult to show that a solution 
to ([55]) exists for every z £ H+. 

The Laplace transform L{z) is equal to 

13'^ if 



,{x,y) lim - — — V {ijj^rjyrjk'nk+i 

rn-oo Zn + 1 ^ — ' 



|fe|<r 



= ^I]E[(/.,(x,2/)((5i(a;-y) + ,5_i(x-y)] 



E 



^(j)^{x,x + 1) 



We define the operator , acting on the set of real valued functions u on , by 

1 



(36) 



{Tzu){x,y) = 



-(Vu)(x,2/) 



Then (|35|) can be written in the following form 

Ipz + Tz4>z = 

where is the function given by 

{5i{x)6o{y) + 5Q{x)5i(y)) 



Pz{x,y) = -- 



Observe that HTj^Hoo < (4/A)||0||oo so that if A > 4 then is contractive for the 
that for A sufficiently large 



norm. It follows 



fc=0 



(37) 



For any x e Z^ , we have the following representation of the operator 

(r^)(x) 

= E ... E (ei-w)...(efc-w)e-^J="^^(^+'=i+-+"^'w(x + ei + 



|eil = l |efcl = l 
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with the convention that the term in the exponential corresponding to j — is 14 (x). We obtain 
(/)^(x,a; + 1) 



— ■^R.w. 



.fe=0 



[6i{x + Sl)So{x + Si) + 5o{x + Sl)Si{x + SD] 



By summing over a; G Z and by taking the expectation with respect to the disorder, we obtain 



L{z) 



1. 



k=0 



fc=0 

By taking first the expectation with respect to 7, we see that we can translate the environment and hence 
the potential by {Sf, Sf) in the first expectation and by {Sl, Sl) in the second one. Therefore, we get 

□ 

Lemma 13. There exists Aq > such that, for any z G _ff+ with ^^(z) > Aq and almost every disorder 7, 
the limit oj Lm{z) as N goes to infinity exists and is given by 

-| 00 

(38) - 2^^"^^'* ['^^^^'^^ '''^'^^S^iiSl - Sl) 

k=0 

Proof. The proof is very similar to the previous one. We look for a solution in the form 

hz iv) =^(l>z{x,y)r]^'ny 
with (l)z{x,y), {x,y) £ T^, the solution of 



x,V 



(39) 



{z + ilx +ly)'i-x^y)<l>z{x,y) + {W (j> z) {x , y) = 



-y) + 5-i{x - y) 



Let A be the real part of z G H+ and define the operator Tz, acting on the real valued functions on T^, 
according to p6p . Then pop can be written in the form 0^ + Tz4>z = Pz, where pz is the function given by 



Pzix,y) = - 



{Si{x-y) + S-i{x-y)) 

2{Z + {j^ + ■^y)la:=iy) 



Observe that ||Tzu||oo < (4/A)||u||oo so that if A > 4 then Tz is contractive for the || • ||oo norm. Therefore 
we have the following representation of (pz 



00 

k=0 
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For any x e we have 
(T^)(x) 



|eil = l 



J2 (ei • w) . . . (efc • w)e- ^?=o ^-(-+<=i+-+^.)m(x + ei 



ek) 



Hence, we obtain 



(r,V.)(x,x + i) = --(-4)* 



Since 14((a;,a;) + Sj) = T^jT^C'S'j), the ergodic theorem implies 



1 1 r 



fe=0 



This completes the proof. 



□ 



Since the sequence {L']^^{z))ff is a bounded sequence of analytical functions on Montel theorem 
implies it forms a compact sequence in the Banach space of analytical functions. Let L'^'^,L'^'^ be any 
(analytical) limit points corresponding to the realizations of 7^ and 7^ of the disorder. For any z G iJ+ such 
that {^{z) > Ao}, we have 

Since the two analytical functions i^o'^ and L^^'^ on coincide on {z; D\{z) > Aq} with , they are 
equal on to . It follows that, for almost every realization of the disorder and every z G the limit 
as N goes to infinity of L^^(z) exists and is equal to L^{z). The theorem is a trivial consequence of the 
following non trivial fact: 

Lemma 14. The limit, as X € (0,+cx3) goes to 0, of L^{X) exists. 

Proof. The proof is similar to the proof of Theorem 1 in [5] (see also [3]). □ 
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